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Abstract 

We present a general scheme to derive higher-order members of the Painleve VI (PVI) 
hierarchy of ODE's as well as their difference analogues. The derivation is based on a 
discrete structure that sits on the background of the PVI equation and that consists of a 
system of partial difference equations on a multidimensional lattice. The connection with 
the isomonodromic Garnier systems is discussed. 
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1 Introduction 



In recent years there has been a growing interest in discrete analogues of the famous 
Painleve equations, i.e. nonhnear nonautonomous ordinary difference equations tending 
to the continuous Painleve equations in a well-defined limit and which are integrable in 
their own right, cf. Even though the qualitative features of the solutions of these 

systems are not yet fully understood, nonetheless in most of the known examples the 
main ingredients of their integrability have been exhibited. Recently, a classification of 
continuous as well as discrete Painleve equations in terms of the root systems associated 
with affine Weyl groups, has been proposed on the basis of the singularities of the rational 
surfaces of their initial conditions and their blowings- up, cf. |p. 

In a recent paper, we established a connection between the continuous Painleve 
VI (Pyi) equation and a non-autonomous ordinary difference equation depending on four 
arbitrary parameters. This novel example of a discrete Painleve equation arises on the 
one hand as the nonlinear addition formula for the Pyi transcendents, in fact what is 
effectively a superposition formula for its Backlund-Schlesinger transforms, on the other 
hand from the similarity reduction on the lattice (cf. @, |5| ) , of a system of partial difference 
equations associated with the lattice KdV family. In subsequent papers, [^, some more 
results on these systems were established, namely the existence of the Miura chain and the 
discovery of a novel Schwarzian PDE generating the entire (Schwarzian) KdV hierarchy of 
nonlinear evolution equations and whose similarity reduction is exactly the Pyi equation, 
this being the first example of an integrable scalar PDE that reduces to full Pyi with 
arbitrary parameters. 

In the present note we extend these results to multi-dimensional systems associated 
with higher-order generalisations of the Pyi equation. Already in |^] we noted that the 
similarity reduction of the lattice KdV system could be generalised in a natural way to 
higher-order differential and difference equations, without, however, clarifying in detail the 
nature of such equations. What we will argue here is that, in fact, such equations constitute 
what one could call the Painleve VI hierarchy and its discrete counterpart. Whilst the idea 
of constructing hierarchies of Painleve equations by exploiting the similarity reductions of 
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hierarchies of nonhnear evolution equations of KdV type is at least two decades old, cf. 
[P, the issue has gained renewed interest in recent years, cf. e.g. [^-[1^], because of the 
hypothetical possibility that these hierarchies of higher-order Painleve equations yield new 
transcendents. Evidence to that effect might be given by the asymptotic analysis of the 
higher-order equations, since they seem to be governed by hyper-elliptic functions rather 



than elliptic ones as is the case for the original Painleve equations, |14|. 

Most of the existing results on hierarchies of discrete and continuous Painleve equations 
are restricted to the examples of Pj and Pn hierarchies, since only in these cases it is clear 
what hierarchies of nonlinear evolution equations should be taken as the starting point for 
their construction. In the case of the other Painleve equations, notably Pyi) it has been 
less clear what to take as a starting point for the construction of its hierarchy. With the 
results of [^, ^, |7| we are now well-equipped to tackle this problem, and in the present paper 
we outline the basic construction of the equations in the discrete as well as continuous 
Pyi hierarchy. In fact, we shall demonstrate that the lattice KdV system can be naturally 
embedded in a multidimensional lattice system achieving the higher-order reductions by 
including more terms in the relevant similarity constraint which provokes the coupling 
between the various lattice directions. 

It should be noted that in a sense higher-order Pyi systems already were constructed 



by R. Garnier in his celebrated paper of 1912, [15|, extending the original approach of R, 



Fuchs who was the first in |lf:] to find Pyi arising from the isomonodromic deformation of 
a second-order linear differential equation. We will conclude our paper with a discussion 
of these Garnier systems, which in view of the recent interest in algebraic solutions of Pyi) 



cf. e.g. [17|-|2C], deserve in our opinion some renewed attention. 



2 The Discrete PVI Hierarchy 

In following earlier work e.g. ^, cf. also |l|], a coherent framework was developed in 
which the similarity reduction of both discrete as well as continuous equations associated 
with the lattice KdV family were treated. Surprisingly, from these reductions the full 
Pyi equation for arbitrary parameters emerged together with a four-parameter discrete 
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equation, i.e. a discrete Painleve equation. From the treatment of Q it was evident how to 
to extend the lattice system of partial difference equations and their similarity constraints 
leading to the reductions to Pyi- Here we describe explicitely this higher-dimensional 
lattice system and discuss their explicit reductions. 

The lattice KdV family of equations contains many related equations such as the 
lattice Schwarzian KdV, the lattice modified KdV (mKdV) and the actual lattice KdV 
equations. We concentrate here on one member of this family only, namely the lattice 
mKdV equation: 

PVn,mVn,m+l + Q Vn,m+lVn+l,m+l = QVn,mVn+l,m + P Vn+l,mVn+l,m+l (2.1) 

cf. e.g. [^], with discrete independent variables n,m and depending on additional pa- 
rameters of the equation p, q, i.e. the lattice parameters. As was pointed out earlier, cf. 
pl| , the lattice equation ( ^.1] ) actually represents a compatible parameter-family of partial 
difference equations: namely, we can embed the equation into a multidimensional 

lattice by imposing a copy of ( |2.1| ) with different parameters on any two-dimensional sub- 
lattice, identifying each lattice direction with a corresponding lattice parameter E C in 
which direction the sites are labelled by discrete variables nj (noting that these are not 
necessarily integers, but shift by units, i.e. Ui £ 0, + 6i G C). Thus, combining two 



different lattice directions, labelled by (i,j) we can write the lattice equation (2.1) on the 
corresponding sublattice as 

Pivv-' + PjV-^v^-' = pjvv^ + piV^v^-' (2-2) 

in which we use the right superscripts i,j to denote the shifts in the corresponding di- 
rections, whereas we will use left subscripts i,j denote shifts in the reverse direction, 
i.e. 

v = v{n;p) , =Tjv{n;p) = v{n + ej;p) , p = Tj'^v{n;p) = v{n - efp) , 

where n denotes the vector of the discrete variables nj, for all lattice directsion labelled 
by i, each corresponding to the component pi of the vector p of lattice parameters. We 
use the vector Sj to denote the vector with single nonzero entry equal to unity in its j^^ 
component. 



The consistency of the lattice equation ( |2.2|) along the multi-dimensional lattice follows 
from the diagram of Figure 1: considering the three-dimensional sublattice with elementary 
directions {61,62,63} then on each elementary cube in this lattice the iteration of initial 
data proceeds along the six faces of this cube, on each of which we have an equation of the 
form (2.2). Thus, starting from initial data v, v^, v'^, we can then uniquely calculate 
the values of w^^, v^^ and v'^^ by using the equation. However, proceeding further there are 
in principle three different ways to calculate the value of w^^^, unless the equation satisfies 
(as is the case for the equation ( |2.2| )) the special property that these three different ways 
of calculating this point actually lead to one and the same value. It is indeed at this point 
that the consistency of the embedding of the lattice MKdV into the multidimensional 
lattice is tested. In fact, by a straightforward calculation we find that this value is given 
by 



{pf - pDpjv'-v'' + (pj - Pi)pkV^v'- + {pI - p'^)piV^v^ 



1 h 3i A; — 1, 2, 3 



{Pi - Pk)Pj'^^ + ip] - p'i)Pkv'' + {pI - p'j)piV^ 

(which is clearly invariant for any permutation of the labels ijk), independent of the way 
in which we calculate this value! Thus, the equation ( p. 2D can be simultaneously imposed 
on functions v{ni,n2,n3, . . . ) of the lattice sites. This is precisely the discrete analogue of 
the hierarchy of commuting higher-order flows of the (modified) KdV equation! 
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Figure 1: Consistency of the lattice equation. 



As a consequence of this compatibility we will call the system (2^) a holonomic system 
of partial difference equations. 



Now, we turn to the issue of the symmetry reduction of the multidimensional lattice in 
the sense of |5|. It follows from the general framework of that the similarity constraint 
for the multidimensional lattice MKdV system is as follows: 

Y,n^ai = ^i-v , zy = A(-1)^''^S (2.3) 

i 

/i and A being constants, and in which the variables are given by 

V'^ — fl} , 

ai = -. . (2.4 

The sum in ( |2.3D is over all the i labelling the lattice directions, the choice of which decides 
the order of the reduction. To analyse the reduction we need a number of relations for the 
objects tti which follow from (^]^), namely 

l+aj = iP^X^^-P^)i-^+^)+^P^ ^ (2.5) 

PiXij + Pj 

^ Pj i-^ij -^ij ~^ Pi{ i-^ij -^ij) Pj 

Pj iXij Xij — Pi{iXij + Xij) + Pj 

Pj i^ij Xij Pj( iXij Xij^ -\- Pj 

Pj iXij Xij ~\~ Pii^ iXij ~t~ Xij 'j -\- Pj 

in terms of the following variables: 



, ^/j (2.6) 



~ yij ' ~ i ~ yj ' ^ ) 

Xij = -J , iXij = T- ^Xij = —J. (2.8) 

The variables and Xij = 1/Xji are not independent, but related via: 

_ PiXij+pj _ -piXij+pj 

j^ij — ■ — Xij — — , (^zJ.yj 

PjXij+pi pjXij-pi 



as well as 



Tj^ Xij 1 dj 



(2.10) 



Xij 1 + ai 

We note that since the left-hand side of (|2.6| ) depends only on the label i but not on j, 
for fixed i this represents a set of — 2 coupled first-order ordinary difference equations 
with respect to the shift in the discrete variable between the A^ — 1 variables Xij, j ^ i. 
Furthermore, the relations ( ^.51 ), for the same fixed label i, provide us with a set of — 1 
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first-order relations between the variables aj, j ^ i, and thus together with the similarity 
constraint (2.3) where Oj is substituted by ( |2.6| ) we obtain a set of 2(A^ — 1) first-order 



nonhnear ordinary difference equations for the 2(A^ — 1) variables Xij, Uj, j ^ i, which 
together form our higher-order discrete system. In the next section we will explicitely 
disentangle this coupled system in the cases N = 2 and = 3. 

The continuous equation for the PVI hierarchy derive from the differential equations 
with respect to the lattice parameters pi, which read: 

d 

-piTT- log V = Uitti . (2-11) 

dpi 



It can be shown that the differential relations ( |2.11| ) are actually compatible not only 
amongst themselves, but also with the the discrete equations on the lattice ( |2.2D , i.e. the 
discrete and continuous flows are commuting: 

d f dv\ d f dv\ dv^ ^ / dv\ 
dpi\dpjj dpj\dpij ' dpj \OpjJ 

This can actually be demonstrated by explicit calculation exploiting the discrete relations 
( |2.5D , but we will not give the details here (which follow closely the pattern of calculations 
of Q). Thus, we have a coherent framework of a large multidimensional system of equa- 
tions with discrete (in terms of the variables rii) as well as continuous (in terms of the 
parameters pi) commuting flows, in terms of which compatible equations of three different 
types (partial difference, differential-difference and partial differential) figure in one and 
the same framework: the partial difference equations are precisely the lattice equations 
( |2.2|) , the differential-difference equations are the relations ( p. 11 ), whilst for the partial 



differential equations in the scheme we refer to our recent paper Q. Here we will focus 



now on the reductions under the symmetry constraint (2^) in order to derive closed- from 
ODE's in terms of the lattice parameter pi. To make this reduction explicit we use ( 2.11j ) 
in combination with (p.5|)-( priO ) to obtain differential relations for the Oj, namely 



da^ = [(1 + a,,)(l - a,)Xj, - (1 + a,)(l - a,)X,,] , (2.12) 

opi Pj - pi 

as well as the following relations for the reduced variables Xij 

+ni ^ 2 {^ji ~ ^ij) + '^fc 2 i-^ik — X^i) 

Pi - Pj Pk - Pi 



in which we have abbreviated 

-Y- _ {piXjj -Pj){pj -piXji) _ ^ 

Pj - pf 

Using ( |2.12| ) in conjunction with ( 2.13| ) and using the similarity constraint to ehm- 
inate the Oj, we obtain a coupled first-order system of differential equations w.r.t. the 
independent variable ti = pf in terms of the 2N — 2 variables Ofc, Xik , {k ^ i). Solving 
the variables from the linear system given by the equations ( |2.13| ) and inserting them 
into we obtain a coupled set of second-order nonlinear differential equations for the 

variables Xi^. 



3 Special Cases: N=2, N=3 

We will now analyse the basic relations of the general framework presented in the previous 
section in the cases = 2, 3 only in order to arrive at slightly more explicit equations, 
demonstrating that the reduction leads to ordinary difference equations (in the discrete 
case) or to ordinary differential equations (in the continuous case). 



N=2: 



We will be very brief about the two-dimensional case N = 2 which was the main subject 
of study in the earlier paper |^ . There the compatibility of the similarity constraint and 



the lattice equation was stated, and the various relations resulting from (2.5) were already 
written down. Using in this case the slightly simpler notation: 

ai = a , a2 = h , xi2 = x , X\2 = X 

and using the discrete independent variables rii = n, n2 = m, as well as the lattice param- 
eters pi = p, p2 = q, we derived the second order nonlinear non- autonomous difference 
equation: 

2(„+l) ^ 2^^^_^^^_jj„^^^^^^ 



1 - yn+lUn 1 - ynUn-l 

^ (^ - A(-l)")(r^ - l)vn + r(l - yl) [{n + \)-{m + |)(-1)"] ^^^^ 

{r + yn){l + ryn) 
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(using a slightly different notation from the one of [^), where r = p/q and where the 
variables y„ are related to the X and x by the prescription: y2n = x{2n) for the even 
sites, whilst y2n+i = —1/ X{2n + \) for the odd lattice sites (the latter choice being mainly 
motivated by the wish to cast the equation into a convenient shape). It was pointed out 
in that whilst a continuum limit of ( p.l| ) yields the Py equation, its general solution 
can be expressed in terms of Pyi transcendents (noting its dependence on four arbitrary 
parameters, /i, A, r and m). 

The continuous equation for the variable X in terms of the lattice parameter p as 
independent variable in this case reads: 

Pip" - q^fX{qX - p)ipX -q)^ = 

= \p{p' - [m(3X2 ^ 1) _ 2(^2 + ^2)^] ' ^ 

+ {q'-p') [2p'X{pX - q){qX-p) + {q^-p'fX^] ^ 

+ \q [{aX^ - (3)ipX - qf{qX-pf + {p" - q^)X^ ((7 - l){qX - pf -{5- l){pX - qf)] , 

(3.2) 

and it is not difficult to show that this is actually the Pyi equation through the identifi- 
cation w{t) = pX{p), where t = p^, and setting (? = 1, leading to 

(fw 1 / 1 ^ 1 ^ 1 \ fdwV /I ^ 1 ^ 1 \ dw 



dt"^ 2 \w w — 1 w — t J \ dt J \t t — 1 w — t J dt 

w(w-l)(w-t)f t-1 ,^ ^t(t-l)\ 

with the identification of the parameters a,P,^,S as follows: 

a = {fx — I' + m — nf' , P = {fj, — u — m + nf' , 

7=(/i + z^ — m — n — 1)^ , 5 = {11 + V + m + n + l)"^ . (3.3b) 

Eq. ( |3.2D is interesting in its own right since it provides us with a covariant way of writing 
Pvi, noting its invariance under the transformations: 

n <-> m , p <-> q , X ^ 1/X . 
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N=3: 



This first higher-order case deals with the first genuinely multidimensional situation of 
three two-dimensional sublattices, on each of which a copy of the lattice MKdV equation 
( |2.2D is defined. In addition there is also the similarity constraint (|2.3| ) which couples the 
three lattice directions. Thus, for the three-dimensional case we have a coupled system of 
equations whose symbolical representation is shown in figure 2. 

In the previous section we have already demonstrated the consistency of the three copies 




Figure 2: symbolic representation of lattice equation and similarity constraint. 

of the lattice equation ( |2.2| ) amongst themselves. What remains is to investigate the 
compatability of the lattice equation and the similarity constraint, demonstrating that the 
determination of the values of the dependent variable v by using the lattice equation in 
all three directions plus the similarity constraint is unique (assuring the single-valuedness 
of the solution around localised configurations). 

In Figure 3 we have indicated how the iteration of the system proceeds starting from 
a given configuration of initial data (located at the vertices indicated by •) and moving 
through the lattice by calculating each point by means of either the lattice equation (points 
indicated by o) or the similarity constraint (points indicated by x). The first point where 
a possible conflict arises, due to the fact that the corresponding values of the dependent 
variable can be calculated in more than one way, is indicated by (8>. It is at such points that 
the consistency of the similarity reduction needs to be verified by explicit computation. 
This has been carried out for this three-dimensional case using MAPLE. Obviously, the 



9 



iteration involves too many steps and the expressions soon become too large to reproduce 
here. 




Figure 3: Consistency of the constrained lattice system. 



In order to analyse the explicit reduction in this case, we redefine the following objects 

ai = a , a2 = b , 03 = c 

X12 = X , = Y , X12 = X , Xi3 = y 

using also ni = n, n2 = m, 723 = h, as well as Pi = p, P2 = Q and P3 = r to simplify the 
notation. To start with the continuous equations, fixing the independent variable to be p 
we obtain the following linear system for the quantities b and c from eq. (|2.13| ) 







mb 




x + y 2y 




he 




(^_^)X + (n + m)^ 





+ + log ^ 

/" + + log y 



+ 



(3.4) 



- + {n + h)^i^ -Y)+ m^(i - X) 
where we have used the similarity constraint to elininate the quantity a. Furthermore 
from ( |2.12| ) we obtain the differential relations 
d{mb) 



dp 

djhc) 
dp 



q 



jYiq 1 

; 7^[{n + fi - - mb - hc)(l - b) — 

' — p'^ X 



-{I + b){n - fi + u + mb + hc)X] 
hr 



1 



^2 _^i[{n + li-y-mb- hc){l - c)- 
— (1 + c){n — ^ + v + mb + hc)Y] . 



(3.5a) 
(3.5b) 
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Solving b and c from the linear system (|3.4| ), and substituting the results in the differ- 
ential relations ( 3.5a| ) and ( 3.5b| ), we obtain two coupled second-order nonlinear ODE's 
depending effectively on six free parameters, namely /i, v, n, m, h and q/r. 

Alternatively, we can derive a system of second-order ordinary difference equations by 
fixing one of the discrete variables, say n = ni, and using the relations (|2.5| ) to obtain the 
equations 



ipX-q)b + pX + q = ' (6 + 1) 



{pY — r)c + pY + r 



qX -p 
{r^-p'^)Y. 



(3.6a) 
(3.6b) 



rY — p 

where the tilde denotes the shift in the lattice direction associated with the variable n. 
Using the similarity constraint 



na + mb + he = II — v , v = X{—1) 



n+m+h 



(3.7) 



to eliminate the variables c, we obtain the following linear system in terms of b and b 

{q^-p^)X -{pX - q){qX - p) 
—m{r'^ — p'^)Y m{pY — r){rY — p) 

(pX + q){qX - p) - {q^ - p^)X 
{rY - p){{pY - r){fi - v - no) + h{pY + r)) - (r^ -p^)Y{h + ^i + v - {n + \)a) 

where the a and a can be expressed in terms of X and Y by 

qX X+p{X -X) - q rYY+p{Y-Y)-r 



(3.8) 



qX X -p{X +X) +q rYY -p{Y +Y) + 



(3.9) 



(where the undertilde denotes the backward shift with respect to the discrete variable). 
The system of equations (3^), (3/7) and (|3.9|) - or, equivalently, (|3.8D together with ( ^ ) 



leads in principle to a fourth order ordinary difference equation in one variable. In fact, 
solving b and b from ( |3.§| ) and then eliminating b altogether by a shift in the independent 
variable n we get a coupled system containing one equation in terms of X, X, X, X and 
Y, and the equation (|3.9| ) which is first order in the both X and Y with respect to the 
shift in the variable n. This system of equations depends effectively on six free parameters, 
namely fi, u, m, h, q/p and r/p. 
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4 Isomonodromic Deformation Problem 



The isomonodromic deformation problem for the multidimensional lattice system is of 



Schlesinger type, |22]. In the two-dimensional case it was already presented in Q for 
special values of the parameters /x, v, cf. also for the general parameter case. The 
extension from the two-dimensional to the multidimensional lattice is immediate: one only 
needs to introduce additional terms of similar form for each additional lattice direction. 
Thus, the Lax representation consists on the one hand of the linear shifts on the lattice of 
the form 

iP'iK) = Tit/jiK) = Li{K)i;{K) , (4.1) 
in which k is a spectral parameter, and where the Lax matrices Li are given by 

Hk) = I , (4.2) 

leading to the Lax equations 

L{Lj = Liu (4.3) 

which lead to a copy of the lattice MKdV equation on each two-dimensional sublattice 
labelled by the indices (^, j). On the other hand we have the linear differential equation 
for iP{k) with respect to its dependence on the spectral variable n 




the compatibility of which with ( [4.1| ) leads to the similarity constraint ( |2.3D . In addition, 
we have differential equations for ijj in terms of its dependence on the lattice parameters 
Pi which are of the form 

dpi Pi \ / + I 1 
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for each of the variable pi. It is the variables U = pf that play the role as independent 
variables in the continuous PVI hierarchy. 

The elimination of the back-shifted vectors T~^ip by using the inverse of the Lax 
relations (|4.1| ) lead to the following linear differential equation for tp 

thus leading to the problem in the Schlesinger form, with regular singularities at 0, oo, {ti}. 
The matrices Ai are given by 



A, 







Ai — 71,; 



The continuous isomonodromic deformation is provided by the linear differential equations 
in terms of the lattice parameters, namely 

where 



n,- 



2k \^ i(l + a,) 

Eq. (^]^) is not quite in standard form, and we need to apply a gauge transformation of 
the form 

ij^VtP , V=\ , (4.8) 

\u/v I I 

to remove the term with Pj, where the auxiliary variable U obeys an interesting set of 
equations by itself (in fact this is the object obeying the lattice KdV system of equa- 
tions) , cf. P], but we will not give any details here. With this gauge, the continuous 



isomonodromic deformation (4/7) adopts the standard form 



^ = , A, = VA,V-' , (4.9) 

whilst the discrete isomonodromic is readily obtained from the Lax representation of ( [4.1| ) . 
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5 Connection with Garnier Systems 



Interestingly, already M.R. Garnier in his seminal paper of 1912, |[T^, embarked on the 
question of finding higher-order analogues of the PVI equation, adopting the method 
that was proposed somewhat earlier by R. Fuchs, in |16], which can be identified with the 



isomonodromic deformation approach, cf. also |22]. Garnier gave a general construction of 
such higher-order equations constituting coupled systems of partial differential equations, 
which are the isomonodromic Garnier systems. As a particular example, he wrote down 



explicitely in |15] the first higher-order PVI equation in terms of the following coupled 
system, consisting of the second order ODE in terms of two dependent variables w = w^t, s) 
and z = z(t, s) 



92 



w 



1 1 
— + — 

w w — 



+ 



1 



1 1 

-'1+— + 



w 

1 



t 



+ 



1 



1 



w 

1 



w 

1 



1 

+- 



1 ' t — s t — w 
w{w — l){w — s){z — t) 



dw 

zl \ dt 

dw 

Zl 'dt 
2 



dz 



w 



+ 



2 z{z -l){z- s){w - t){z -w) \dt J {z-t){z- w) 
2w{w - l){w - t){w - s){z - tf 



dw\ / dz 
'dtj [di 



i2(t-l)2(t 



s)2(u) 



a + f3 + -f + 6 + K + 
t{t - l){t - s) 



7 tsa + l (t_i)(s_i) p+l 



+ 



z [z 
s{s - l)(s - t) 



(z-i) {w-iy 



{z -t) {w- ty 



+ 



{z — s) {w — sy 



(5.1a) 



together with coupled first order PDE's 

t{t — 1) dw ^ s{s — 1) dw 



w{w 



t- z dt 



s — z ds 



w 



t{t-l)dz s{s-l)dz z{z-l) 



+ 



z — w 



(5.1b) 
(5.1c) 



t — w dt s — w ds 

It should be pointed out that the system consisting of ( ^.la ), ( ^.Ib ) and ( 5.1c| ) amounts 
actually to a fourth order ODE in terms of w = w{t) only, and as such can be rightly 
considered to be the first higher-order member of the Painleve VI hierarchy. In fact, 
Garnier gave in his paper a number of important assertions: i) that his system of equations 



14 



is completely integrableQ, and that it degenerates (under the autonomous limit) to a system 
living on the Jacobian of a hypereliptic curve, ii) that the symmetric combinations of the 
dependent variables of the system, as functions of each one of the essential singularities 
(i.e. singling out one of the independent variables) are meromorphic in terms this variable 
except for the fixed critical points which are at 0,1, cxd, or at the location of the values 
of the other independent variables^, in) that for the parameters of the system in general 
position the symmetric functions of the dependent variables are essentially transcendental 
functions of the constants of integration (i.e. of the initial data). 

Subsequent work on the Garnier systems was done mostly by K. Okamoto and his 
school, cf. e.g. |2^, However, it seems that in most of these works these systems were 
treated rather as an underdetermined system of PDE's rather than (as Garnier himself 
clearly had in mind) as a consist system of ODE's. Although it is not easy to find the 
explicit transformation of the lattice system exposed in sections 2 and 3 to the systems 
that Garnier wrote down, in particular to find the explicit relation between the above 



system ( |5.1D and the system consisting of (3^) and (p^), it is to be expected that such a 
mapping exist. The identification is probably easiest to obtain via the transformation of 
the corresponding Schlesinger type of system as given in section 4 and the linear system 
that Garnier exploited in |15]. However, the search for such an identification will be left 
to a future study. 

Let us finish with some remarks on the relevance of these results for work that is been 
done in recent years. One of the most exciting developments is the way in which the issue 
of algebraic solutions of Pyi have arisen in recent years, e.g. in connection with WDVV 



equations, Frobenius manifolds and quantum cohomology, cf. e.g. the review [25]. Such 
algebraic solutions were already known to Picard, Painleve and Chazy. In fact, in his early 
paper |lf:] R. Fuchs obtained a realisation of Pyi in terms of an elliptic integral, and this 



realisation was subsequently used by Painleve in [26|, to derive an elliptic form for the 

^Obviously, Garnier's use of the term integrability was meant here in the precise sense of that of a 

compatible system, very much in the same sense as the compatibility of the continuous and discrete 

systems that we have encountered in sections 2 and 3. 

^This assertion amounts to the well-known Painleve property. 
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Pvi equation, a form of the equation that was recently recovered by Manin in |. The 



assertions of Garnier in [15| on his generahsation of the Fuchs' approach might form a 
starting point for extending this ehiptic connection to the Garnier systems, in which case 
we would expect to be able to find a realisation of those systems in terms of hyperelliptic 
integrals rather than elliptic ones. This might eventually lead to the construction of 
algebraic solutions of those systems, possibly in the spirit of the recent papers [|^, 28 1. It 
would be of interest to further investigate the role of the discrete systems in connection 
with the Garnier systems: we expect them to constitute the superposition formulae for the 
underlying higher root systems of the corresponding affine Weyl groups. Thus, eventually, 
a geometric interpretation of the Garnier systems and their discrete analogues in the sense 
of the blowings-up of the corresponding rational surfaces of their initial conditions, along 
the lines of the recent paper M, might be anticipated. 
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